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Abstract. 

A case of long-range transport of a biomass burning plume from Alaska 
to Europe is analyzed using a Lagrangian approach. This plume was sam- 
pled several times in the free troposphere over North America, the North At- 
lantic and Europe by 3 different aircraft during the IGAC Lagrangian 2K4 
experiment which was part of the ICARTT /ITOP measurement intensive 
in summer 2004. Measurements in the plume showed enhanced values of CO, 
VOCs and NO y , mainly in form of PAN. Observed O 3 levels increased by 
17 ppbv over 5 days. A photochemical trajectory model, CiTTyCAT, is used 
to examine processes responsible for the chemical evolution of the plume. The 
model was initialized with upwind data, and compared with downwind mea- 
surements. The influence of high aerosol loading on photolysis rates in the 
plume is investigated using in-situ aerosol measurements in the plume and 
lidar retrievals of optical depth as input into a photolysis code (Fast-J), run 
in the model. Significant impacts on photochemistry are found with a de- 
crease of 18% in O 3 production and 24% in O 3 destruction over 5 days when 
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including aerosols. The plume is found to be chemically active with large O3 
increases attributed primarily to PAN decomposition during descent of the 
plume towards Europe. The predicted O3 changes are very dependent on tem- 
perature changes during transport, and also, on water vapor levels in the lower 
troposphere which can lead to O3 destruction. Simulation of mixing/dilution 
was necessary to reproduce observed pollutants level in the plume. Mixing 
was simulated using background concentrations from measurements in air 
masses in close proximity to the plume, and mixing timescales (averaging 
6.25 days) were derived from CO changes. Observed and simulated O3/CO 
correlations in the plume are also compared in order to evaluate the photo- 
chemistry in the model. Observed slopes changed from negative to positive 
over 5 days. This change, which can be attributed largely to photochemistry, 
is well reproduced by multiple model runs even if slope values are slightly 
underestimated suggesting small underestimation of photochemical processes. 
The possible impact of this biomass burning plume on O3 levels in the Eu- 
rope boundary layer is also examined by running the model for a further 5 
days, and comparing with data collected at surface sites, such as Jungfrau- 
joch, which showed small O3 increases and elevated CO levels. The model 
predicts significant changes in O3 over this 10 days due to photochemistry 
but the signal is largely lost due to the effects of dilution. However, measure- 
ment in several others BB plumes over Europe show that O3 impact of Alaskan 
fires can be punctually significant over Europe. 
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1. Introduction 

There is a large body of evidence that pollutants can be transported from source regions 
to remote areas, and thus can influence the oxidizing capacity over downwind regions, and 
occasionally regional air quality thousands kilometers away from emissions. During export 
from emission regions, polluted plumes undergo chemical evolution and mixing with air 
masses of different origin so that their chemical composition can radically change. Ozone 
levels in such polluted plumes are of particular interest since 0 3 is a major pollutant, 
damaging for health but also the major source of the OH radical, and the most efficient 
cleansing agent in the atmosphere. O 3 is mainly produced by reactions involving oxidation 
of volatile organic compounds (VOCs), and carbon monoxide (CO) in the presence of 
nitrogen oxides (N 0 iT =N 0 +N 02 ), and mainly destroyed by reactions involving H 2 O and 
HO s radicals. 

O 3 evolution in anthropogenic pollution plumes transported over long distances has 
been widely studied [Wild and Akimoto, 2001; Jacob et al, 1999] but less is known about 
0 3 formation in biomass burning plumes, especially those in boreal regions. Biomass 
burning (BB) has been recognized as an important source of many trace gases and parti- 
cles [Crutzen et al., 1979; Crutzen and Andreae, 1990] including CO, C0 2 , CH 4 , VOCs, 
oxygenated organic compounds and NO z . BB emissions occur most extensively in tropical 
and subtropical regions with the largest contributions from Africa and South America, 
and smaller contribution from southeast Asia, Indonesia and Australia [ Galenter et al, 
2000]. Boreal forest fires in North America and Russia are less frequent, but impor- 
tant in summer time in certain years [Stohl et al, 2006]. Long-range transport of BB 
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plumes has been previously observed during several campaigns. For example, African 
BB plumes were sampled over the tropical south Atlantic [Andreae et al, 1994] and over 
Australia [Pak et al, 2003], southeast Asian BB plumes were sampled over the west Pa- 
cific [Rondo et al, 2004] whilst Canadian BB plumes have been sampled over Tennessee 
(United States) [Wotawa and Trainer, 2000; McKeen et al, 2002] and also at the surface 
station Mace Head (Ireland) [ Forster et al, 2001]. All these plumes showed enhanced CO 
correlated with enhanced O 3 . Modeling studies of BB plumes (either with box models or 
global models) also often show net O 3 production in such plumes [Mauzerall et al, 1998; 
Duncan et al, 2003; Tang et al, 2003]. 

The IGAC Lagrangian 2K4 experiment was conducted in summer 2004 over the North 
Atlantic as part of ICARTT [ Fehsenfeld et al, 2006], also encompassing the European 
Intercontinental Transport of Ozone and Precursors (ITOP) experiment. The aim of this 
experiment was to perform a Lagrangian study of polluted air masses, making intercep- 
tions several times using 4 different research aircraft during their transport from North 
America to Europe. Lagrangian experiments have been conducted previously, for example, 
in the stratosphere [Rex et al, 1998] using balloons measurements, or close to emissions 
regions [Daum et al, 2000; Volz-Thomas and Kolahgar, 2000]. This is the first time that 
a Lagrangian experiment has been conducted in the free troposphere on intercontinental 
scales. The main advantage of this kind of study is that the exact chemical change in a 
plume over several days can be estimated by comparing two Lagrangian samplings, and 
the processes leading to these changes can be evaluated with reduced uncertainty in the 
initial chemical composition of the plume. After the campaign a detailed analysis us- 
ing trajectories, a Lagrangian dispersion model, and in-situ measurements was performed 
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by Methven et al. [2006], and showed that several Lagrangian samplings of the same air 
masses had indeed been achieved. The period of the campaign was characterized by ex- 
tensive fires in Alaska (the largest on record), and northern Canada [ Damoah et al. , 2006; 
Stohl et al. , 2006]. The quantity of CO released during this period is believed to be of the 
same order as anthropogenic emissions from the entire continental United States for the 
same period [Pfister et al., 2005; Turquety et al, 2006]. CO enhancements clue to these 
fires were observed together with aerosols, and trace gas enhancement in the Arctic [Stohl 
et al, 2006], in North America, both in the boundary layer (BL) and the free troposphere 
(FT) over Chebogne Point (Nova Scotia) [ Duck et al, 2006], and Texas [Morris et al, 
2006], in the lower free troposphere over the Atlantic Ocean [ValMartin et al, 2006], 
and in the FT over Europe [Petzold et al, 2006; Ravetta et al, 2006]. Deep convection 
associated with these fires (so-called pyroconvection) may have also led to injection of 
BB plumes into the upper troposphere (UT), and even into the lower stratosphere (LS) 
[Damoah et al, 2006]. 

In the study reported here, we focus on one particular case of the transport of a BB 
plume over the North Atlantic that was sampled by 3 different aircraft over a period 
of 5 days (section 2). To evaluate the relative contributions of chemical and dynamical 
processes to the evolution of O 3 concentrations in the plume during transport, a photo- 
chemical trajectory model has been used (section 3). BB plumes are also characterized by 
high aerosol loadings, and therefore, the model was improved to include the effects of BB 
aerosols on photolysis rates by making use of measured parameters to estimate aerosol 
optical properties (section 4). This is the first time, to our knowledge, that measured 
parameters have been used to constrain photolysis rate calculations over intercontinen- 
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tal distances. In section 5, results are presented from model runs initialized with upwind 
data, and compared to downwind measurements. The model was run along one trajectory 
representative of the plume transport pathway, and which gave the match between plume 
samplings. The model was run first with chemistry only, and then including mixing with 
air masses in close proximity to the plume. The ability of the model to reproduce the 
contribution of photochemistry to O3 changes in the plume, was further tested in section 
6, by examining the evolution of O3/CO correlations. Model runs initialized with mea- 
surements taken across the plume sampling were compared to downwind 03:C0 ratios. 
Finally, we examine other available data collected downwind over Europe in order to see if 
a signal of this BB event could be detected. Further model runs were used to estimate the 
overall change in O3 concentrations since the upwind sampling 10 days earlier (section 7). 
Whilst the results presented here are for one particular event, such detailed Lagrangian 
analyzes are very valuable since they allow direct evaluation of processes occurring in a 
plume during long-range transport. 

2. Lagrangian Forest Fire Plume Case 

Methven et al. [2006] used a novel technique to identify Lagrangian matches between 
flight segments from different aircraft during the entire IGAC Lagrangian 2K4 campaign. 
This technique combined Lagrangian trajectories from the United Kingdom Global Atmo- 
spheric Modeling Project (UGAMP) model using global meteorological ERA-40 analyzes, 
hydrocarbon fingerprint analysis, and results from a particle dispersion model (FLEX- 
PART). A match was defined as an occasion when a pair of whole air samples collected 
during different flights exhibited highly correlated hydrocarbon fingerprints (when such 
measurements were available), and the sample-time windows were also connected by both 
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backward and forward trajectories. Results from FLEXPART run with CO tracers were 
also used to confirm matches. Five clear Lagrangian cases covering a variety of situations, 
such as low and upper level export of pollutants from North America, and long-range 
transport of BB plumes from Alaska and Canada over the North Atlantic to Europe 
were identified. An important consideration is whether instrument errors can lead to 
observed differences between samplings. Intercomparison flights between the different 
aircraft showed good agreements for O3, and differences of a few ppbv for CO. Further 
details, as well as a critical evaluation of this Lagrangian technique, can be found in 
Methven et al. [2006]. 

In this paper, we focus on one of these Lagrangian cases, the case of long-range transport 
of an intense BB plume sampled over 5 days first by the National Aeronautics and Space 
Administration (NASA) DC8 over Newfoundland on 18 July, then by the U.K. Facility 
for Airborne Atmospheric Measurements (FAAM) BAel46 flying out of the Azores on 20 
July, and lastly the German Deutsches Zentrum fr Luft- und Raumfahrt (DLR) Falcon 
flying west of Paris on 23 July. The British and German aircraft were flying as part of 
the European ITOP project. The segments of the three flights identified as Lagrangian 
matches were labeled as Lagrangian case 2 by Methven et al. [2006]. 

Backward trajectories from the match segments of the 3 flights clearly show that the 
sampled plume had an origin in northern Canada and Alaska (see Fig. 1) where forest 
fires were burning during this period [ Pfister et al., 2005]. The FLEXPART model, run 
with daily CO forest fire emissions from Canada/ Alaska, showed that BB plumes ob- 
served during the campaign were usually an accumulation of biomass burning emissions 
over several days which were then transported into the free troposphere either by frontal 
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uplift or convection [Stohl et al. , 2006]. In this case, FLEXPART results show that the 
pollutants sampled by the aircraft were primarily emitted by fires on the 14 and 15 July 
(see FLEXPART products at http : // www.al.noaa. gov /icartt .analysis/). 

In this study, a photochemical trajectory model has been used to analyze observed O3 
changes. Therefore, it was necessary to select a trajectory representative of the plume 
pathway between the upwind DC8 segment (used to initialize the runs), passing via the 
BAel46 match, and ending at the Falcon match. Trajectories from the FLEXTRA model 
[Stohl et al, 1995] were used here, but it should noted that differences between trajectories 
from this model, and those calculated by Methven et al. [2006] are small since they used 
the same wind fields (European Centre for Medium-range Weather Forecast (ECMWF) 
ERA 40 analyzes), and for this particular case, the models produce the same features. 

As noted previously, Methven et al. [2006] used a combination of forward, and back- 
ward trajectories, trace gas and VOC analysis, and FLEXPART results to match up 
plume segments. Use of trajectories alone shows some differences. For example, forward 
trajectories from the DC8 flight segment arrive around 1km lower than the Falcon match 
segment. Specific humidity, taken from ECMWF, is also a factor of 2 wetter. Conversely, 
backward trajectories from the Falcon segment are about 1km above the DC8 match, and 
the analysis temperature slightly colder. Errors of this nature are inherent in trajectory 
calculations [Stohl, 1998] gave positioning errors on the order of 20% of the travel dis- 
tances). Runs along both backward and forward trajectories, and also examination of 
the sensitivity of O3 change to temperature and water vapor (see Section 5.2) led us to 
use a backward trajectory from the Falcon plume match segment which gave the best 
agreement with observed meteorological variables. 
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The chemical evolution of the plume can be investigated first by examining the differ- 
ences in the chemical data averaged over the time match windows (see Table 1 and Table 
1 in Methven et al. [2006]). Pollutant levels measured by the DC8 were representative 
of the chemical composition of the BB plume 3 to 4 days after emission. By this time, 
CO values measured by the DC8, at a mean altitude of 7.5 km, were still exceptionally 
high with a mean value of 449 ppbv, and maximum values up to 600 ppbv (see Fig. 2) 
which correlated with high values of several hydrocarbons such as ethene, ethane and 
acetylene. NO,, was mainly in the form of PAN with relatively low levels of HNO3 and 
NO x observed in the plume. Unfortunately PAN was only measured at the beginning of 
the plume sampling (plume edge) by the DC8, and so is not representative of the mean 
plume value which are likely to have been higher in the plume center. 

O3 levels in the plume were rather low (mean value of 62 ppbv) compared to values 
outside the plume. However, Pfister et al. [2006] simulated background O3 levels lower 
than 40 ppbv over the Alaskan source region during summer 2004 with a global model. 
Very low background (non-fire affected) values (20 to 30 ppbv) have also been measured 
previously in this region during ABLE-3A, so that the mean value of 62 ppbv observed 4 
days later by the DC8 may indicate that O3 production has in fact taken place although in- 
situ measurements in the emission region would be needed to confirm this. Interestingly, 
O3 and CO were anti-correlated in the DC8 plume (see section 6 for further discussion). 

By the time the BAel46 sampled the plume 2 days later, CO had decreased by about 
30 ppbv, and by the time it was sampled by the DLR Falcon southwest of Paris, 3 days 
later, CO had decreased to half the levels seen by the DC8. In contrast, O3 levels had 
increased by 17 ppbv over 5 days between the different samplings. In this study, we aim 
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to explain this O 3 increase. It is also interesting to examine whether a plume like this 
with elevated O 3 can have an impact on O 3 levels in the lower troposphere or boundary 
layer over Europe. This is investigated in section 7. 

3. Modeling tools 

The Cambridge Tropospheric Trajectory model of Chemistry And Transport (CiTTy- 
CAT) was used to examine the different processes influencing the evolution of trace gases 
within the plume during transport, and in particular the evolution of O 3 and CO. CiT- 
TyCAT was initialized with measured values in the upwind plume on 18 July, and run 
along the trajectory representative of the plume pathway. Results from the model were 
then compared with downwind measurements taken in the matching flight segments of the 
BAel46 and Falcon flights, 2 and 5 days later. In order to take into account the influence 
of BB aerosols on photolysis rates, the CiTTyCAT model was coupled with the Fast-J 
photolysis scheme [ Wild et al . , 2000] which can calculate photolysis rates in the presence 
of aerosols. 

3.1. Photochemistry trajectory model: CiTTyCAT 

CiTTyCAT has been previously used to examine the origin of polluted layers over the 
North Atlantic during previous campaigns [Wild et al., 1996; Evans et al., 2000]. The 
model is considered as an isolated air parcel, and run along trajectories calculated using 
large-scale meteorological analyzes. Eighty-eight chemical species are treated in the model 
including degradation of 12 hydrocarbons using chemical rate data from JPL (2003), and 
updates discussed in Arnold et al. [2006]. The photolysis scheme normally used in the 
model is based on a 2-stream multiple scattering scheme [Hough, 1988]. Emissions and 
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dry deposition can be included in the boundary layer as well as wet deposition in presence 
of clouds. A scheme for mixing with other air masses can also be included. As discussed 
previously, the model was run along a trajectory calculated by the FLEXTRA model 
using ECMWF windhelds. Temperature and water vapor used in the chemical scheme 
were also interpolated from the large-scale ECMWF analyzes along the trajectory. In this 
study, the possible effects of heterogeneous reactions on BB aerosol have not be taken 
into account. At present, the impacts of such reactions on this kind of aerosol, consisting 
primarily of organics, is very uncertain. Whalley et al. [2006], who also examined this 
case, attempted to take this into account by including the loss of N2O5 on BB aerosol 
with a 7 of 0.02 as suggested by work from Evans and Jacob [2005]. They found a 20% 
reduction in NCR concentrations over the Erst 24 hours which resulted in negligible O3 
changes compared to runs without heterogeneous loss. Wet deposition processes are not 
considered here as aerosol measurements on board the Falcon as well as ECMWF analyzed 
cloud fields suggest that the plume did not pass through any clouds. 

Table 1 gives the mean plus/minus the standard deviation of all the measurements 
collected by the aircraft used in the model initialization (DC8 data - except PAN) and 
comparison (BAel46 and DLR Falcon data). The total number of measurements used 
to calculate these values in the match segments is also given. As explained in section 2, 
measurements were missing during the main part of the DC8 match so the measured 
PAN / CO correlation during the BAel46 match was used to estimate initial PAN values. 
This correlation was scaled by measured DC8 CO as also discussed in Whalley et al. [2006] 
giving a mean PAN value of 2800 pptv together with a 400 pptv standard deviation. This 
assumes that PAN and CO do not change much chemically between the two matches. 
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Indeed, as will be seen in section 5.2, simulated CO and PAN only decrease, due to 
chemistry, by 20 ppbv and 200 pptv, respectively, producing only a small impact on the 
PAN/CO correlation. 

3.2. Photolysis scheme: Fast-J 

The Fast-J scheme calculates photolysis rates (j-values) in the presence of prescribed 
cloud and aerosol layers [Wild et al, 2000] making it possible to simulate directly the 
physical properties of scattering and absorbing particles including the full, untruncated 
scattering phase function and the total optical depth. For the case under investigation, 
aerosol optical properties (AOPs) were derived from Falcon measurements in the BB 
plume. BB optical depths (ODs) used in the code were calculated from a combination of 
DC8 and Falcon retrievals for the 18 and 23 July as discussed in the next section. The 
original 2-stream photolysis scheme in CiTTyCAT uses 171 wavelengths intervals in the 
short-wave whereas Fast-J solves the 8-stream multiple scattering problem by using only 
7 wavelength bins selected in the range 289-850 nm to minimize errors in key j-values over 
a range of atmospheric conditions. Values of jNC >2 and j 03 (i.e. that producing 0 1 D) 
calculated using the two codes run in CiTTyCAT with no clouds or aerosols, were compa- 
rable to within 1%. In runs with no aerosols or clouds discussed in the following sections, 
a simple profile of absorbing black carbon uniformly distributed from the surface to 10 
km with an optical depth of 0.001 was included to represent relatively clean conditions. 

4. Aerosol measurements and optical properties 

BB aerosols consist mainly of black carbon (BC), organic carbon (OC) and potassium. 
It is the BC component which acts as a strong light absorbent. Sampling obtained during 
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the campaign in a BB layer on board the BAel46 showed that these aerosols were princi- 
pally made of organic and black carbon with the organic fraction dominating (Coe, pers. 
comm.). In the plume considered here, analysis by Petzold et al. [2006] suggests a BC 
fraction of 3.8% by volume compared to typical BC mass fractions in temperate/boreal 
fires on the order of 2 to 7% [Reid et al, 2004], 

Previous studies have shown that BB aerosols can significantly decrease photolysis rates, 
and thereby influence 0 3 evolution in BB plumes [Tang et al, 2003; Trentmann et al, 2003; 
Duncan et al, 2003]. However, their impact depends strongly on AOPs and aerosol con- 
centration. BB aerosol AOPs can vary significantly depending on the age of the aerosols, 
the mode of combustion of the fire (smoldering or flaming), the combustion efficiency, 
and the combustible material [Reid et al, 2005]. Optical properties of BB aerosols which 
are wavelength dependent or their absorption properties can be very different from one 
fire to another. For example, aerosols from temperate and boreal fires often have higher 
values of the single scattering albedo than aerosols from tropical fires. Indeed, a realistic 
evaluation of the influence of BB aerosols on photolysis rates requires a good evaluation 
of the AOPs from particular fire events. 

In-situ measurements of BB aerosol characteristics have been used previously to simulate 
aerosol impacts on photolysis rates, and thus on photochemistry, in the initial stages of 
plume transport during the first few hours since a fire [Trentmann et al, 2002], To our 
knowledge, this is the first time that a combination of in-situ, and also airborne /ground- 
based lidar measurements have been used to estimate the time evolution of AOPs over 
several days, and their subsequent impact on photolysis rates and photochemistry during 
long-range transport. 
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4.1. Calculation of Aerosol Optical Properties 

During the campaign, the DLR Falcon aircraft was equipped with extensive in-situ 
aerosol instrumentation. The data has been used to derive particle size distributions, 
particle number densities, and particle absorption coefficients. More information on the 
aerosol measurements can be found in Petzold et al. [2006] . Aerosol measurements in the 
BB plume on the 23 July were used to derive the refractive index, and size distribution of 
the BB aerosols following the procedure described in Fiebig et al. [2002], This information 
was then used as input in a Mie code to derive BB AOPs [Mishchenko et al., 1999]. 

Firstly, a composite size distribution in the plume was derived from the measurements, 
and a three-modal logarithmic normal distribution was fitted to the particle size distribu- 
tion (see Figure 9 in Petzold et al. [2006]). Each of these three modes were parameterized 
by integral particle number concentration, N, particle median radius, r m , and the geomet- 
ric standard deviation, sigma (<r 9 ), as given in Table 2. 

For the refractive index calculation, an internal mixture of ammonium sulfate and soot 
was used with a soot fraction of 3.8%. As noted previously, a large fraction of the aerosol 
likely consisted of organic carbon, but here ammonium sulfate is used as a surrogate 
since it has similar refractive index [ Sloane , 1984], The soot fraction of 3.8% was derived 
from the measurements [ Petzold et al, 2006] so that the calculated absorption coefficient 
matched the measurements as closely as possible and the hypothesis of internal mixture 
was also confirmed by the measurements. Finally, the calculated refractive index using 
these hypotheses was of 1.54 +0.0168L 

With these assumptions about the aerosol size distributions, and the refractive index, 
the extinction efficiency (Q), the single scattering albedo (ir 0 ), and the first 8 terms of 
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the Legendre expansion of the scattering phase function (■ w l , i=0,7) were calculated for 3 
wavelengths characteristic of ultraviolet (UV) and visible wavelengths (300, 400 and 600 
nm) using a Mie code. These results were directly used in the Fast-J code to simulate BB 
aerosol optical properties and are summarized in Table 3. In the selected BB plume, w o, 
which represents the fraction of scattering versus light absorption by aerosols, is in the 
upper range of previous estimates derived from measurements in aged boreal and North 
American fire plumes (0.85 to 0.91 over wavelengths from 400 to 700 nm) [Reid et al, 
2005], 

4.2. Calculation of Aerosol Optical Depth 

AOPs alone are not sufficient to describe the effect of a BB plume on photolysis rates. 
Information about aerosol concentrations, and their vertical distribution are also required. 
This information is given by the optical depth (OD) of the aerosol layer which corresponds 
to the integration of the extinction coefficient over a layer. During long-range transport 
of an aerosol layer, the OD changes with time as the layer becomes dispersed. Here, two 
ODs were derived from in-situ measurements of the DC8 aerosol lidar on the 18 July and 
from DLR Falcon on the 23 July [Petzold et al . , 2006]. On the 18 July an OD of 0.2 
at 1024 nm and 0.3 to 0.4 at 550 nm was evaluated for the BB layer with an estimated 
vertical dimension of 7 km which corresponded to a strong increase in the backscatter 
signal from the DC8 lidar (see Fig. 3). For the 23 July, a layer thickness of 2 km was 
used to calculate the aerosol optical depth at 550 nm using an average aerosol absorption 
coefficient of 6.5 x 10~ 6 m -1 , and a single-scattering albedo of 0.9 at 550 nm. The height 
of the layer was taken from vertical profile data collected by the Falcon spanning 3.5 to 
5.5 km. The aerosol absorption coefficient was determined from PSAP (Particle Soot 
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Absorption Photometer) data assuming that the layer-averaged value of the absorption 
coefficient corresponds to 75% of the peak value inside this layer which was 9 x 10 -6 m -1 
[Petzold et al, 2006]. These values yield an aerosol optical depth of 0.13 for the BB layer 
probed by the Falcon on 23 July. 

In addition, a dual wavelength depolarization aerosol lidar operated by the SIRTA 
Observatory in Palaiseau [ Haeffelin et al, 2005], southwest of Paris at about the same 
location as the Falcon measurements, reported approximately 1 hour before the Falcon 
sampling, a layer extending from 2.5 to 5.5 km with high aerosol loading. The lidar data 
was analyzed using a multi-test algorithm applied to all the lidar channels. The layer 
consisted of a lower layer between 2.5 and 4.0 km characteristic of BB aerosols (1064/532- 
nm color ratio of about 0.25 and 2% depolarization), and an upper layer between 4.0 and 
5.5 km containing larger particles (0.5 color ratio). No measurements were available later 
on this day, but it is likely that this BB layer was part of the same BB layer observed by 
the Falcon. An optical depth of 0.1 to 0.15 was estimated for the entire layer based on the 
lidar power attenuation between the observed lidar return above the aerosol layer, and 
the theoretical molecular lidar return in the absence of aerosols [Cadet et al, 2005]. This 
confirms that, even if the Falcon may not have sampled the entire vertical distribution 
of the BB plume, the Falcon profile was representative of the main vertical extent of this 
aerosol layer. 

4.3. Photolysis scheme evaluation 

In order to evaluate the performance of the Fast-J code for this case, simulated and 
measured jN02 values were compared by running the code along the DC8 flight track on 
the 18 July (see Fig. 2). This flight was chosen because a clear signal was seen in the 
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measured data when the aircraft passed through the BB aerosol layer (as seen by the on- 
board aerosol lidar), and measured photolysis rates were reduced. In a run with no clouds 
or BB aerosols the model significantly overestimates the data between 18hl5 and 21h30. 
This time window corresponds to an increase in the measured upward jNC>2 flux and high 
aerosol scattering ratios observed by DC8 backscatter lidar at low levels indicating the 
presence of low cloud. The agreement is partially improved by including a low cloud layer 
(see Fig. 2, blue rectangle) with an OD of 5 extending from 0 to 0.5 km between 18hl5 to 
21h30 (continuous line). However, measured jNC>2 is still overestimated in the BB plume 
(see coincident CO data). A further run was performed including a BB aerosol layer with 
an OD of 0.3 (see Fig. 2, red rectangle) co-located with the high aerosol scattering ratio 
signal measured by the DC8 lidar between 4 and 11 km (Fig. 3). AOPs were prescribed 
as discussed in the previous sections. The 30% reduction in jN02 observed in the BB 
plume (dashed line) is now simulated well by the model suggesting that the AOPs used 
in the Fast-J code are reasonably realistic. Observed photolysis rates of other species like 
O3 or PAN are also reduced in the BB plume, and well simulated by the run including 
the effects of BB aerosols. The impact of these changes on O3 levels are discussed in the 
next section. 

5. Results 

In this section, results are discussed from CiTTyCAT simulations initialized with data 
on 18 July and run for 5 days. Firstly, the impact of high aerosol loading on photolysis 
rates, and photochemistry in the plume are discussed. Then, the influence of both chem- 
ical processes, and dynamical processes related to mixing with other air masses during 
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transport, are examined, and results compared to observations collected downwind during 
the Lagrangian matches. 

5.1. Aerosol impact on photolysis rates 

Concerning the inclusion of BB aerosols in Fast-J for this case, several assumptions, 
and simplifications had to be made. The aerosol layer was assumed uniform in aerosol 
concentration, and its thickness was varied from 7 km on the 18 July to 3 km on 23 July 
with the OD decreasing linearly from 0.3 to 0.13 over 5 days. It was also assumed that 
the trajectory was situated in the middle of the BB plume. The DC8 lidar data (see 
Fig. 3) confirms that the DC8 match is, indeed, in the middle of the BB plume, and the 
CO vertical profile measured by the Falcon on the 23 also shows that the Falcon match 
is approximately located in the middle of the plume. 

In this study, BB aerosols were considered as an internal mixture of black carbon, 
an absorbent component, and sulfate, which has scattering properties. In an absorbing 
aerosol layer, available light is reduced, and therefore photolysis rates decrease whereas 
they can be enhanced in a scattering aerosol layer. For a mixture, the relative influence of 
scattering versus absorption depends on the solar zenith angle, SZA. Scattering properties 
outweigh absorption properties at low SZA, but their influence decreases with SZA. This 
dependence of the scattering/absorption with SZA is dominated by the shape of the 
aerosol phase function which can vary strongly depending on the refractive index, and 
the size distribution of the aerosols. SZA dependence was strong for the BB aerosols in 
this case, and at lower SZA the photolysis rates were even enhanced when aerosols were 
included. Consequently, measured photolysis rates should be analyzed carefully because 
changes due to the presence of BB aerosols will depend on SZA, and therefore on the 
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time of the day when the measurements were made. Another important aspect of the BB 
AOPs is their dependence on the size parameter 2nr/X. For small particles, this leads 
to an inverse dependence of the single scattering albedo on wavelength in the UV. This 
means that for a single particle, smaller wavelengths will be absorbed less. However, a 
plume contains thousands of particles which interact with each other leading to more 
complex effects which depend, not only on the single scattering albedo, but also on SZA, 
the optical depth, and the aerosol phase function [ Reuder and Schwander, 1999]. 

The BB aerosol impact obtained by averaging over the 5 day simulation was found to 
be a reduction in all photolysis rates, and in particular, a reduction of -21% for jNC >2 
compared to -22% for j0 3 . This can be compared with results from Duncan et al. [2003] 
for Indonesian fires where monthly averaged reductions in jNC >2 and j0 3 in the FT over 
source regions were around -10 to -20% using a mean AOD of 0.3. Simulated photolysis 
rate reductions in a Mexican fire transported into Arizona were larger (-20 to -60%), 
but this study used a lower value of the single scattering albedo (0.82) than used here 
[. Kotamarthi et al, 2006]. Also, in contrast to our simulation, jN0 2 was reduced more 
than j0 3 because different size distribution and refractive index were used. This shows 
that results are very sensitive to the AOPs parameters used for particular fire events. 

In this case, the reduction in photolysis rates leads to a slowing down of the photo- 
chemistry with a reduction of -18% in 0 3 production and -24% in 0 3 destruction. 0 3 
production is reduced mainly because lower jN0 2 leads to reductions in NO plus H0 2 . 
0 3 destruction is mainly due to the reaction of O x D with water vapor, and this term 
is reduced the most when aerosols are included due to the reduction in j0 3 . But, this 
term, and also the second most important 0 3 loss term (reaction of 0 3 with HO 2 ) are also 
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reduced due to lower O3 levels resulting from less O3 production. The overall effect of 
including these effects on net 0 3 production depends on the balance between production 
and destruction. Changes in photolysis rates due to the inclusion of BB aerosols will re- 
duce net O3 production in a plume mainly producing O3, and increase net O3 production 
in a plume mainly destroying O3. It will be seen in section 7 that after 8 days of simula- 
tion, the plume enters a net O3 destruction regime. From this moment, the reduction in 
photolysis rates leads to an increase in 0 3 production. 

5.2. Chemistry only simulation including aerosol radiative effects 

Simulated concentrations over 5 days of CO, 0 3 , NOy, NO, PAN and water vapor are 
shown in Fig. 4, and compared to the range of measured values along the matching 
Lagrangian segments on the 20 July (BAel46), and the 23 July (DLR Falcon). In order 
to better characterize the variation of the chemical composition in the plume, 3 chemical 
initializations were used. One using mean observed DC8 values (black lines), and two 
others using mean ± standard deviation (std) values (red and blue dashed lines). It 
should be noticed that the red (blue) simulation corresponds to initialisation with the 
mean +std (-std) for species like CO, VOC, PAN but the mean -std (+std) for 0 3 , since 
0 3 is anticorrelated with precursors such as CO. 

The first important discrepancy in the results is found in the simulated and observed 
concentrations of CO. The small CO decrease observed during the first 2 days is well 
reproduced by the model but after 5 days CO is overestimated by more than 150 ppbv. 
NO y is also overestimated by about 1000 pptv after 5 days. As will be seen in the 
next section, these overestimations are likely to be due to dilution of the plume during 
transport. 
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Despite the large discrepancies in CO and NO,,, it is still interesting to examine this 
chemistry-only simulation in order to understand the photochemical processes in such BB 
plumes during long-range transport. The large O3 increase, of about 17 ppbv observed over 
5 days between the DC8 and the Falcon matches, is reproduced, and even overestimated 
by the chemsitry-only simulation giving a mean O3 increase of 22 ppbv in 5 days just due 
to photochemistry. 

Diagnosis of model results for these chemistry-only simulations shows that the large 
increase in O3 is mainly due to the high PAN loading in the plume. PAN is stable at low 
temperatures encountered during the first 2 days of the simulation, but rapidly decomposes 
at high temperatures encountered once the plume subsided leading to release of N0 2 , and 
resulting in increased O3 production. Unfortunately downwind PAN measurements are 
not available but modeled NO values fall within the range of the observations on both 
on the 20 and 23 July. Note that in these runs, NO,, levels in the plume are maintained 
because part of the N0 2 released from PAN decomposition is transformed into HNO3 
which makes up about 30% of the NO,, level on 23 July. 

In order to illustrate the effect of temperature on PAN decomposition, several stationary 
steady state runs at a fixed position corresponding to the location where the DC8 sampled 
the plume, were performed. Mean DC8 measurements, together with scaled PAN, in the 
plume were used for the chemical initialization, and only the temperature varied between 
each run. Values of PAN and O3 after 1 day of simulation are shown on Fig. 5. At 
temperatures lower than 260 K no changes in PAN are seen. PAN thermal decomposition 
begins at temperatures higher than 260 K, and is non-linear with temperature especially 
beyond 270K. At 290K, 90% of the initial PAN is decomposed after 1 day. Therefore, 
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increasing temperature higher than 290 K will not have a strong effect as most of the 
PAN has already been decomposed thus limiting the potential for further O3 production. 
It is interesting to note that the value of 90% PAN decomposition after 1 day at 290K is 
not dependent on the initial PAN value. This also explains why, in the simulation, even if 
the simulated temperature increases by 25 K (from 240 to 265 K) between the DC8 and 
the BAel46 matches, only a small fraction of PAN is decomposed. 

In model runs with and without PAN, it was found that PAN decomposition accounts 
for 80% of the total O3 increase between DC8 and Falcon matches. In the run without 
PAN, conversion of the initial NO loading into O3 through NO2 photolysis is only sufficient 
to maintain net O3 production during the first 2 days resulting in net O3 destruction over 5 
days. Since PAN has a central role in governing O3 production, it is important to consider 
uncertainties in the PAN decomposition rate which can vary by ± 100% at 280K accord- 
ing to the latest IUPAC recommendations ( http : / /www.iupac — kinetic.ch.cam.ac.uk ). 
Increasing this rate by a factor of 2 leads to faster PAN decomposition, lower PAN (by 
500 pptv), and higher O3 (by 6 ppbv) concentrations on 23 July. However, this represents 
a 10 ppbv overestimation of the Falcon data. 

It is clear that physical variables such as temperature and water vapor are also im- 
portant in governing the production, and destruction of O3 during long-range transport 
of such plumes. Comparison of observed values in the plume, and values interpolated 
from ECMWF analyzes are shown in Fig. 4. Analyzed and measured water vapor, as 
given by specific humidity agree well. Nevertheless, it is interesting to see how O3 levels 
respond to changes in water vapor. Water vapor in the plume is low during the 5 first 
days so that halving the water vapor content in the model has a very small impact on 
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O 3 net production (less than 1 ppbv). On the contrary doubling the water vapor content 
along the trajectory increases 0 3 destruction, and leads to a decrease of 18% in net 0 3 
production. 

Measured and simulated temperatures agree reasonably well except at the beginning of 
the simulation where temperatures are underestimated by more than 7K on the 18 July, in 
part because the trajectory used passes slightly above the DC 8 match. However, increasing 
the initial temperature to match the measured one does not change 0 3 production as both 
temperatures are too low for PAN decomposition. Measured temperatures during the 
BAe-146 and the Falcon matches are similar but ECMWF analysis show that between 
these two matches, the temperature reached 278K (at the southern most point along 
trajectory), which leads to the large 0 3 production. A sensitivity test with a constant 
temperature between these two matches gives a decrease of 41% in net 0 3 production, 
and therefore much lower 0 3 after 5 days (9 ppbv less). 

These sensitivity studies show how important specific humidity and temperature are 
for 0 3 chemistry. Therefore, it is important that meteorological parameters used by 
models agree as closely as possible with observations. The evolution of meteorological 
parameters along the plume transport pathways is also crucial for estimating 0 3 levels, so 
that photochemical evolution cannot be inferred from measured values alone. It will also 
be shown (in section 7) that water vapor and temperature play a key role in determining 
0 3 levels in plumes transported into the lower troposphere where high water vapor can 
lead to strong 0 3 destruction balanced by 0 3 production through PAN decomposition. 

Overall, whilst chemical changes alone could explain observed 0 3 changes during plume 
transport, they cannot explain the observed changes in CO and NO,,. 
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5.3. Simulation of mixing/dilution of the BB plume 

In the previous section, model runs with no mixing were discussed, and the plume was 
assumed to be isolated from the background and not diluted. In the real atmosphere, a 
plume is subject to stirring by large-scale winds followed by mixing. This mixing process 
entrains surrounding air masses, and thus changes concentrations in the plume [Stohl, 
2001]. The scheme for mixing in CiTTyCAT is a simple linear relaxation scheme [ Evans 
et al., 2000] resulting in an exponential decay of plume concentrations towards back- 
ground concentrations. The rate of mixing is defined by the e-folding timescale K rmx in 
day -1 which can also be written in terms of a mixing timescale r in day (r = jJ - — ). 
This scheme requires the estimation of so-called background concentrations as well as the 
mixing timescale. In this study, instead of trying to estimate average background concen- 
trations, concentrations were defined using in-situ measurements in air masses adjacent 
to the BB plume. The mixing timescale was chosen in order to obtain the best agreement 
between measured and modeled CO which is a good tracer for mixing since it is relatively 
long-lived on the timescales considered here, and not strongly affected by chemistry. A 
strong decrease in CO is thus characteristic of mixing/ dilution of the plume. 

Analysis of observed CO in the plume showed a small CO decrease between the DC8 
and the BAel46 matches, and a much stronger one between the BAel46 and the Falcon, 
when the plume reached western Europe (decrease of 40% of BAel46 CO). Trajectory 
analysis showed that just south of the BAe-146 match, trajectories reached a stagnation 
point and switched from descent in quiescent clear air to being near a warm conveyor belt 
moving up from the southwest that was likcley to be more turbulent. This resulted in a 
change in direction in the transport pathway (see Fig. 1). It is thus realistic to simulate 
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two mixing regimes representative of these different situations: a slow mixing timescale 
(20 days) during the two first days of simulation following by a faster one (5 days), leading 
to an integrated mixing timescale of 6.25 days. 

The values of these mixing timescales, chosen to match the CO measurements, are 
dependent on whether the same parts of the plume were sampled since, in reality, strong 
gradients exist at plume extremities. Lagrangian matches are not perfect and it is possible, 
for example, that the Falcon may have missed the center of the plume sampled by the 
DC 8 explaining the large difference in measured CO. However, other parts of the same 
large plume were sampled over Europe by different platforms and the measured CO was 
in the same range as the Falcon data or lower (see section 7). It is also worth noting that, 
if the same mixing rates are applied to aerosols in the plume (i.e. mixing the DC 8 OD 
with a zero background), this leads to a decrease from 0.3 to 0.128 after 5 days, which 
is consistent with the estimated OD on the 23 July from in situ measurements and lidar 
data. This supports the choice of the mixing rates used in this study. 

Concerning background values, Methven et al. [2006] showed that this BB plume was 
situated within a dry intrusion that had descended from close to the tropopause over North 
America. DC 8 measurements of CO, O 3 and data from the O 3 and backscatter lidars also 
clearly show that the BB plume was in close proximity to air with low CO and high O 3 
concentrations, characteristic of the UT. It is therefore reasonable to imagine that the 
BB plume and UT air masses have mixed. Therefore, initially, the DC 8 measurements in 
the UT air mass were used to constrain background concentrations (see Table 4) together 
with a mixing timescale of 20 days during the first 2 days and 5 days during the last 3 


DRAFT 


October 26, 2006, 6:25pm 


DRAFT 



REAL ET AL.: PHOTOCHEMISTRY IN ALASKAN FOREST FIRE PLUME 


X- 27 


days. In results from this run (not shown) O3 was overestimated by more than 7 ppbv 
compared to the Falcon measurements. 

Examination of ECMWF water vapor suggests that there was the mixing of wetter 
air into the plume where the trajectories change direction (see water vapor on Fig 4). 
This corresponds to when the BB plume lost contact with drier air above and to north, 
and came into contact with more moist air from the south associated with a warm con- 
veyor belt. Evidence for this split between the BB plume and UT air masses is also 
seen in satellite water vapor images for the 20 and 21 July (Cooper, pers. comm., see 
http : / 1 www. al.noaa.gov /csd/metproducts/icartt200A/). Therefore, a second run was 
performed where both the mixing time scale and background concentrations were changed 
after 2 days in an attempt to represent changes in the composition of air masses mixing 
into the BB plume. In fact, air masses sampled by the Falcon next to the plume on the 
23 July were typical of free tropospheric (FT) air masses, and had higher water vapor 
content and lower O3 than the UT air mass. These measurements were used to constrain 
background concentrations during the last 3 days of the simulation (see Table 4). 

Results using UT/FT background concentrations are shown in Fig. 6. Simulated O3 
and CO agree very well with Falcon measurements giving some confidence in the hypoth- 
esis that the BB plume mixed first with a UT air mass before mixing with a FT air mass. 
As both UT and FT mixed air are poor in PAN and HNO3, NO y in the plume is reduced 
by about 1000 pptv leading to very good agreement with the Falcon measurements. The 
final breakdown in modeled NO,, is 55% PAN, 3% NO x and 42% HNO3. Simulated NO 
values are still in the range of measurements because NO level are governed by photo- 
chemistry in the plume (PAN decomposition), and are not strongly influenced by mixing. 
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It is interesting to note that mixing not only changes concentrations in the plume by 
introducing background air, but can also change the photochemistry in the plume. In this 
case, the main effect is dilution of PAN in the plume leading to a reduction of net O3 
production in the plume by 10%. 

Overall, these results show that inclusion of mixing processes, albeit in a simple manner, 
produce results in better agreement with the full range of measurements. 

6. O3/CO correlations 

In the previous section, it was shown that a combination of photochemistry and mixing 
with UT /FT air masses can explain the evolution of the BB plume since the results agree 
reasonably well with the observations for several trace gases. However, it is very difficult 
to know the real composition of mixed air masses, and thus to evaluate the relative 
contributions of photochemistry and mixing. For example, mixing of the plume with a 
UT air mass could lead to a large O3 increase in the plume even if the plume does not 
undergo any photochemical production. It is therefore possible that both photochemistry 
and mixing simulations could be incorrect but that the final results could match the 
observations because the effects of the different processes compensate each other. However, 
the proposed scenario can be tested in more rigorously by looking at the evolution of trace 
gas correlations, and more specifically the slope of such correlations. Here, we focus on the 
O3/CO correlation since O3 can change rapidly over several days due to photochemical 
processes, and CO is a good tracer of mixing over this time period. Indeed, if mixing in the 
plume is uniform, then O3/CO correlations are not changed by mixing processes and, in 
the absence of loss processes, any change in the slope can be attributed to photochemical 
processes alone. This assumes that all parts of the plume are mixed with the same 
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background, and with the same mixing timescale at any point in time as was assumed 
here in Section 5.3. Considering that the plume evolution is followed during only 5 
days, and that the segment matches correspond to inside the plume and not the edges, 
this hypothesis of uniform mixing seems reasonable. Mixing processes could also have 
indirect effect on O3/CO slopes (hOs/dCO) by diluting O3 precursors, and so modify 
photochemical processes in the plume. 

The case studied here is particularly interesting to look at in terms of O3/CO correla- 
tions since a 4 day old BB plume was measured by the DC8 with a negative dO^/SCO 
value, an unusual feature in an ageing BB plume (see Fig. 7). Negative values of hOs/5CO 
have been measured in young BB plumes due to titration [Stith et al., 1981; Hobbs et a/., 
2003] but most BB plumes observed during previous campaigns showed rapid 0 3 produc- 
tion, and positive <503/<5C0 values after just 1 or 2 hours [Goode et al ., 2000; Andreae 
et al. , 1988]. SO^/SCO measured by the BAel46 in the BB plume was still negative, but 
with a less steep slope than the DC8 (i.e. 6 days after emission). Aged BB plumes with 
negative correlations have also been observed at the PICO-NARE station on the Azores 
during the ICARTT campaign ValMartin et al. [2006]. The reason for these negative 
correlations in aged BB plumes is not clear. It does not necessarily mean that the plume 
has not undergone any photochemical O3 production but that the O3 production was not 
strong enough to change the initial negative slope to positive. High aerosol loadings over 
the source regions [Stohl et al., 2006; Duck et al., 2006] might have strongly dampened 
initial O3 production in these plumes or else other processes such as O3 loss by reac- 
tion with organic component present in aerosols or nighttime chemistry including the of 
heterogeneous hydrolysis of N 2 0 5 [ValMartin et al., 2006] might have played a role. 
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In contrast, the dO^/SCO measured by the Falcon was positive (see Table 5) with a 
value of 0.078. This change from negative to positive slope after 5 days indicates that 
the plume is strongly chemically active as already suggested by the model results dis- 
cussed previously. The slope value measured by the Falcon is in the lower range of values 
reported for previous campaigns which range enormously from 0.045 to 0.88 [ Andreae 
et al. , 1988, 1994; Goode et al. , 2000]. It is also lower than values reported by ValMartin 
et al. [2006] for BB plumes with positive correlations observed during the ICARTT cam- 
paign at the PICO-NARE station ranging from 0.14 to 0.89. However, direct comparison 
of SO3/SCO is difficult because, in our study, the slope calculation is made in the La- 
grangian sampling of the plume which does not include, for example, the edges of the 
plume, and represents only CO values larger than 150 ppbv (for the Falcon match). Us- 
ing PICO-NARE measurements, Pfister et al. [2006] showed that if CO thresholds used 
to calculate <503/<5C0 slopes in a plume are increased from 100 to 120 ppbv, then the 
slope value is almost divided by a factor 2 (from 0.44 to 0.26). This shows that 5O3/5CO 
values are very sensitive to the range of values selected in the calculation. 

Therefore, in order to test the scenario put forward in the previous section, and in 
particular, the photochemistry in the model, multiple model simulations were conducted, 
initialized with DC8 measurements averaged every 30s during the match segment. Mod- 
eled and measured O3/CO correlations and slopes were then compared downwind. With 
this method, each run represents a different part of the plume (black crosses in Fig. 8) 
observed by the DC8, and thus the gradient across the plume. Simulation and compari- 
son of 5O3/5CO are relatively common with global models but, to our knowledge, this is 
the first time that such simulations have been performed with a photochemical trajectory 
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model. Results from the multiple simulations after 2 (red triangles) and 5 (blue diamonds) 
are shown for chemistry only runs (Fig. 8 a), and chemistry plus mixing (UT/FT) runs 
(Fig. 8 b). Observed O 3 /CO correlations are represented by linear fits to the data (con- 
tinuous lines). Measured and simulated slopes on the 20 and 23 July are also summarized 
in Table 5. Overall, the model simulations reproduce very well the change of JO 3 /JCO 
from negative for the DC 8 and the BAel46 samplings to positive for the Falcon match. 
Values of the O 3 /CO slope are relatively well simulated but slightly underestimated. The 
results show that the main effect of mixing is to translate simulated O 3 /CO correlation 
lines towards background CO and O 3 levels leading to better agreement with the range 
of measured concentrations. Only small changes are found between O 3 /CO slopes for 
simulations with (slope of 5.1%) and without (slope of 5.2%) mixing. The slope is slightly 
decreased in the mixing run due to dilution of O 3 precursors which decreases O 3 pro- 
duction. In the following discussion all calculated slopes are from chemistry plus mixing 
simulations. 

There are two explanations for the small discrepancy between measured and simulated 
JO 3 /JCO: the 3 matches are not perfect, and the limits of the matching plume are slightly 
different, or the model slightly underestimates the production of O 3 per quantity of CO, 
even if it reproduces the main features of the photochemical evolution of the plume. 
Concerning the possibility that different parts of the plume were sampled, it can be seen 
in Fig. 8 b that the measured CO values during the Falcon match are more dispersed than 
those simulated by the model. This can signify that the limits of the Falcon match are too 
large. For example, reducing the Falcon match window by 2 minutes leads to a reduction 
in the measured slope from 7.8% to 6.1%, which is closer to the simulated value. 
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Concerning the hypothesis of an underestimation in the modeled O3 production, two 
hypotheses have been examined: an increase in photolysis rates due to the presence of 
clouds, and an increase in the PAN decomposition rate (see Section 5.2). Satellite images 
show that low clouds were present over Europe from 22 to 23 July, and Falcon jN02 data 
collected on these two days often showed a strong variability with some very high values, 
typical of broken cloud fields. To see if an increase in the J-values due to the presence 
of low clouds could increase «50 3 / AGO, sensitivity tests using different cloud ODs on 22 
and 23 July were performed. JO3/JCO after 5 days of simulation were compared (see 
Table 5 for details). 5 O3/JCO increases slightly when increasing cloud OD, up to 6.3% 
with a cloud OD of 20. An increase by a factor two in the PAN decomposition rate leads 
to a final slope of 6.5%. In both cases, simulated JO3/JCO are closer to the measured 
one but mean simulated 0 3 levels in the plume are now overestimated by about 3 to 6 
ppbv. Only a change in the mixing scenario could lead to better agreement with the mean 
plume value but this would require mixing with air masses with much lower 0 3 and no 
evidence for such air masses was found in the measurements made near the plume. In the 
case of doubling the PAN decomposition rate, NO y values in the plume are also slightly 
underestimated. It is finally difficult to conclude on the reason for the small difference 
between simulated and observed slopes. It is likely due to combining effect of non perfect 
match window and small underestimation of photochemistry due to either increased in 
photolysis rates over clouds (which are not simulated) or underestimation in the PAN 
thermal decomposition rate. 

Overall, examination of changes in 0 3 /C0 slopes has allowed more rigorous validation of 
the contribution of photochemistry to O3 changes in this BB plume since the model is able 
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to reproduce the change from negative to positive slope. These results give confidence in 
the modeled photochemical contribution, and support the results presented in Section 5.3. 
They confirm that significant O3 production took place during the plume transport, and 
that mixing had a large influence on CO concentrations. 

7. Possible influence on O 3 levels in the European lower troposphere 

The Lagrangian analyses presented in the previous sections have shown that, for the 
particular case under investigation there was a strong enhancement in O3 due to photo- 
chemistry after 5 days. Other studies have examined the influence of boreal forest fires 
during this period on the O3 budget over the North Atlantic. For example, Pfister et al. 
[2006] estimated a 10% contribution to the tropospheric burden at the end of July 2004. 
Cooke et al. [2006] also estimated an O3 increase of 6 ppbv in the boundary layer over 
the United Kingdom from 21 to 23 July 2004 due to BB emissions, but did not find 
any measurements to support their findings. However, global models often underesti- 
mate concentrations of pollutants in plumes due lack of horizontal and vertical resolution 
[Bousserez et al., 2006; Pfister et al., 2006], and may therefore underestimate the impact 
of BB events downwind over Europe. In this section, we extend the previous analysis 
further downwind in order to investigate the possible influence of the case examined here 
on O3 in the lower troposphere, and the boundary layer over Europe. In the future, it 
would be interesting to examine other cases, and to compare these results with global 
models. 

For the period 23 to 26 July we have examined data collected at surface sites to see 
whether there is evidence of transport of the BB plume studied here into the lower tro- 
posphere or boundary layer over Europe. FLEXPART simulations, run with CO BB 


DRAFT 


October 26, 2006, 6:25pm 


DRAFT 



X - 34 


REAL ET AL.: PHOTOCHEMISTRY IN ALASKAN FOREST FIRE PLUME 


tracers, showed that a large part of the plume was transported northwards over Scan- 
dinavia with a significant fraction transported also over central Europe, and finer scale 
filaments transported southwards over France and the Mediterranean. Data was exam- 
ined from Observatoire de Haute Provence (OHP) (43°55N, 5°42E) in southern France 
where two O 3 lidars and a backscatter lidar were operated during the 1TOP campaign 
(see Ravetta et al. [2006] for details); Jungfraujoch (lat:46°55N, 7°98E) at 3.6 km altitude 
in the Swiss Alps where measurements of many trace constituents including CO, O 3 and 
several hydrocarbons were made as well as meteorological variables, and aircraft profiles of 
CO and O 3 from the Measurement of OZone, water vapor, carbon monoxide and nitrogen 
oxides by In-service Airbus airCraft (MOZAIC) programme [Nedelec et al, 2005]. 

Analysis using back trajectories identified the following features which appear to be 
linked to the Lagrangian case discussed in the previous sections (see Fig. 9). Even if the 
different samplings are not directly linked in the Lagrangian sense, they were probably 
initially part of the same large plume. 

A MOZAIC flight over Paris on the 23 July confirmed the presence of the layer seen by 
the DLR Falcon on the same day. Ozone lidar data collected at OHP on the 24 July showed 
a layer at 3 to 4km which was linked to the larger plume, if not exactly the same part of the 
plume sampled earlier by the Falcon. Measurements made during the period 24 July to 26 
July at Jungfraujoch are shown in Figure 10 together with results from the FLEXPART 
BB tracer simulations when observed CO concentrations exceeded 200 ppbv. Trajectories 
analysis show that on the 24 July these enhancements were linked to a part of the larger 
plume when it was transported into the lower troposphere over Europe. Trajectories also 
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suggest that this air mass continued to descend into the European boundary layer over 
eastern Europe and Mediterranean. 

In order to evaluate the possible contribution of this BB plume on O3 levels over Europe 
the runs discussed in section 5 were run for a further 5 days from 24 to 28 July, and 
compared to the data from OHP, Jungfraujoch and MOZAIC (Paris). The results are 
shown in Figure 11 over 10 days. For the last 5 days of simulation, the aerosol OD 
was fixed at 0.13 and a mixing rate of 5 days. FT background concentrations were 
used to constrain background concentrations as the plume was mainly transported in the 
FT, and no measurements were made in the proximity of the plume after the Falcon 
match. Simulated CO agrees reasonably well with the observations at Jungfraujoch when 
mixing is taken into account (pink shaded area) even if O3 concentrations at Jungfraujoch 
are slightly overestimated. It is also interesting to note that the modeled O3/CO slope 
continues to increase up to 11% until the plume enters a destruction regime when slope 
starts to decrease. Observed and modeled 0 3 /C0 slopes at Jungfraujoch on the 24 July 
are very similar, 7.5% and 7.7%, respectively. Agreement with observed O3 at OHP, where 
O3 levels were higher (around 80 ppbv), is less good supporting the idea that a different 
part of the plume was observed over southern France. The proximity of a tropospheric 
fold above the layer (also seen by the OHP lidar) suggests that additional mixing may 
have taken place with air masses of upper tropospheric origin. 

Bearing in mind the caveats that this is not a truly Lagrangian analysis, and that 
mixing and aerosols conditions were less well constrained for the last 5 days of the plume 
evolution, it is still possible to use the results to examine the relative contributions of 
photochemistry and mixing to O3 levels in the plume after 10 days of transport. As 
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noted above, results from the run including mixing give better agreement with the data 
further downwind, especially over Jungfraujoch, and suggest that the net change in 0 3 
after 10 days was very small or even negative. However, mixing processes mask the impact 
of photochemistry considerably by diluting the signature of the BB plume over Europe. 
Results from a chemistry-only simulation (blue shaded area) show that the plume was 
in a net 0 3 production regime during the first 8 days of simulation leading to a increase 
of 30 to 40 ppbv over this period, primarily as a result of PAN decomposition. During 
the last 2 days there is net O3 destruction over Europe when PAN is almost completely 
decomposed due to elevated temperatures and water vapor has increased strongly. The 
net change in 0 3 for the chemistry only simulation over the 10 day period is 8 to 12 pbv. 

Despite the dilution effect on pollutant levels in the plume, it is clear that BB plumes 
can have an impact on CO and O3 levels in the European lower troposphere and even 
the boundary layer since the data and model simulations show high levels of both species 
over Europe between the 22 and 26 July. Several other BB plumes were also sampled over 
Europe during the ITOP campaign. For example, the Falcon sampled BB plumes several 
times with enhanced CO (up to 350 ppbv), 0 3 (up to 100 ppbv) and NO y (up to 2500 
pptv) including one low level case on 26 July off the coast of France at 1.5 km when CO 
was 270 ppbv and O3 was 95 ppbv. MOZAIC profiles collected during the campaign period 
also showed many plumes with enhanced CO (up to 300 ppbv) and 0 3 (up to 85 ppbv) 
over Paris, Frankfurt and Vienna between 3 and 7 km which are likely to be BB plumes 
with such high CO concentrations. All these plumes give different O3/CO slopes ranging 
from 4% to 40%. This likely reflects their different histories in terms of temperature and 
humidity which, as shown previously, can have a large impact on net photochemical O3 
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production. This is particularly true for such BB plumes which have high PAN loadings. 
These results show that it is necessary to be very careful when attributing O 3 changes 
to photochemistry alone since mixing with other air masses can also have a significant 
impact. It may also makes detection of O 3 anomalies linked to long-range transport of 
pollutants into the European lower troposphere or even the boundary layer extremely 
difficult. 

8. Conclusions 

A case of a BB plume transported from Alaska to Europe, has been analyzed using a 
Lagrangian approach. This plume was sampled three times by different aircraft during 
the IGAC Lagrangian 2K4 campaign (summer 2004) over the northeast North America, 
the North Atlantic and Europe. Lagrangian match segments along the flight tracks were 
defined by Methven et al. [2006] using a combined trajectory/ VO C fingerprint analysis. 
Measurements in the plume show an O 3 increase of 17 ppbv over 5 days together with 
a significant decrease in CO, and a change from negative to positive O 3 /CO correlations 
over the same time period. 

The CiTTyCAT model has been used to study the processes governing the chemical 
evolution of the plume, and in particular, changes in O 3 . The model was initialized with 
upwind measurements in the BB plume, and compared with downwind measurements. 
The influence of high aerosol loadings in the plume on photolysis rates during long-range 
transport was taken into account by combining in-situ aerosol measurements in the BB 
plume and lidar measurements with a photolysis module (Fast-J), implemented in CiTTy- 
CAT. The impact of BB aerosol on photolysis rates was -21% for jN 02 and -22% for j 03 
with the overall impact leading to a reduction of photochemistry in the plume: O 3 pro- 
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duction was reduced by 18% and 0 3 destruction by 24% leading to less net 0 3 production 
overall. Factors such as single scattering albedo, optical depth and phase function were 
found to be important. These results show that such BB aerosol effects on photochemistry 
can be important and should be taken into account in chemical model simulations. 

Despite the slowing down of photochemistry in the plume due to the aerosols, the model 
predicts that the plume was very chemically active. Net 0 3 production was estimated to 
be 22 ppbv over 5 days with a chemistry-only run. This large 0 3 production in the plume 
was driven by thermal decomposition of PAN which was very elevated in the plume, and 
made up 80% of NO y during the upwind sampling by the DC8. Temperatures, which 
were low initially, warmed by more than 30K along the plume pathway leading to this 
PAN decomposition. Low water vapor content (measured) in the plume also helped to 
maintain this high net 0 3 production. 

It was also shown that the chemistry-only simulation was not able to reproduce observed 
values of certain pollutants, especially CO and NO,, after 5 days, and that mixing/dilution 
of the plume must also have been important. In-situ measurements collected in close 
proximity to the plume were used to constrain background concentrations, and mixing 
timescales were chosen to obtain the best agreement between modeled and observed CO 
in the plume in relation to the meteorological situation of the flow. The best results were 
obtained by using concentrations representative of a UT air mass during the first 2 days, 
and FT concentrations during the last 3 days with mixing timescales of 20 days and 5 
days, respectively (giving an integrated mixing timescale of 6.25 days). With this mixing 
scenario, simulated CO, 03, NOy and NO levels were in the range of observations made 
in the plume after 2 and 5 days. 
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Results using multiple model simulations were also used to test more rigorously the 
contribution of photochemistry to observed changes by comparing simulated and observed 
O 3 /CO slopes in the plume. The change from negative to positive correlation during 
transport across the North Atlantic was well simulated by the model even if the observed 
slopes were slightly underestimated. These results confirm that the large increase in O 3 
was due to photochemistry. Extension of the model runs for a further 5 days over Europe 
showed that once the plume entered the lower troposphere, 0 3 destruction dominated. 
Dilution effects were also increasingly important, leading to only very small simulated O 3 
changes after 10 days. These results are supported by comparison with available ground- 
based and aircraft data which showed evidence for this plume being transported into the 
lower troposphere over Europe. 

Overall, the influence of Alaskan BB plumes on O 3 levels over Europe may have been 
particularly important at least in the FT as shown by the different BB plume samplings 
from aircrafts and ground based stations and as suggested by global model simulations 
[Pfister et al, 2006; Cooke et al, 2006]. However, 0 3 changes in such plumes can vary 
a lot since they depend on its photochemical and physical history which will depend a 
lot on PAN loading, temperature and water vapor evolution along its pathway as well 
as dilution effects. The results presented here show that it is important to distinguish 
between these different processes. 

Recent studies predict that boreal wildfires may increase in the future due to climate 
change [ Stocks , 1998], and possibly increased direct human impact [ Mollicone et al, 2006]. 
As forest fire plumes are likely to have high PAN loading due to the abundance of oxy- 
genated volatile organic coumpounds (o-VOC) [Singh et al, 1994; Fan et al, 1994], their 
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ability to transport O3 on intercontinental scales is large, and their impact is likely to be 
important over the entire Northern Hemisphere. 
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DC8 Lagrangian 

BAel46 Lagrangian 

Falcon Lagrangian 


sampling (24 mn) 

sampling (5 mn) 

sampling (14 mn) 

(ppbv) 

440 (68) [129] 

414 (53) [300] 

244 (35) [840] 

0 3 (ppbv) 

61.6 (2.6) [133] 

67.3 (1.2) [300] 

79 (3.6) [840] 

NO (pptv) 

115 (65) [105] 

49 (16) [45] 

29 (6) [840] 

N0 2 (pptv) 

28 (10) [105] 



PAN (pptv) 

1842* (140) [30] 

2727 (197) [300] 


HN0 3 (pptv) 

84 (10) [87] 



NOj, (pptv) 



1956 (243) [840] 

C 2 H 4 (pptv) 

2156 (360) [20] 

1123 (371) [131] 


C 2 H 6 (pptv) 

3595 (475) [20] 

4030 (371) [131] 


C 2 H 2 (pptv) 

1365 (169) [20] 

1123 (132) [131] 


C 3 H 8 (pptv) 

790 (132) [20] 



C 3 H 6 (pptv) 

84 (20) [20] 



C 4 H 8 (pptv) 

12 (3) [20] 



C 4 HiO (pptv) 

220 (40) [20] 

82 (9) [131] 


C 5 H 8 (pptv) 

4 (0.8) [20] 



C 5 H!2 (pptv) 

70 (14) [20] 



C 6 H 6 (pptv) 

710 (128) [20] 

858 (64) [131] 


C 7 H 8 (pptv) 

51 (15) [20] 

133 (17) [131] 


F Tjbmp (K) 

247 (0-2) Mbs 

r 2$, 6 m] 

273 (2) [131] 

H 2 0 (g/kg) 

0.24 (0.04) [133] 

0.29 (0.01) [300] 

1.9 (0.4) [131] 


Table 1. Measurements taken during Lagrangian match windows defined in Table 1 


of Methven et al. [2006]. Data for each window is shown as “mean (standard deviation) 
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N (part /cm 3 ) 

r m (micrometer) 

°9 

mode 1 

1400 

0.037 

1.6 

mode 2 

1100 

0.125 

1.3 

mode 3 

15 

0.175 

1.8 


Table 2. Parameters of the 3 log-normal distribution functions fitted to the DLR falcon 
aerosols measurements. N is the integral particle number concentration, r is the particle 
median radius and a is the geometric standard deviation. 


Table 

of 1.54 


A,nm 

Q 

w 0 

w° 

w 1 

w 2 

w 3 

w 4 

w 5 

w 6 

w 7 

300 

2.6700 

0.9160 

1 

2.143 

2.528 

2.236 

1.688 

1.127 

0.686 

0.369 

400 

1.9490 

0.9140 

1 

2.014 

2.040 

1.480 

0.892 

0.447 

0.232 

0.114 

600 

0.8600 

0.8950 

1 

1.689 

1.330 

0.671 

0.331 

0.162 

0.097 

0.063 

3. Aerosol optical proper 

ties at 3 wave] 

lengths (calcul 

lated with a refractive index 


+ 0.0168i at 550 nm) 



CO 

o 3 

NO, 

NO 

UT/LS 

94 

105 

900 

120 

FT 

105 

70 

550 

8 


Table 4. Concentrations measured in air masses in close proximity to the BB plume 
and used as background concentrations for the mixing simulation (see text for details) 
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Measured 

without 

with 

with cloud 

cloud 

cloud 

PAN rate x 2 



mixing/ clouds 

mixing 

(OD=5) 

OD=10 

OD=20 


Slope / BAel46 

-0.00875 

-0.0151 

-0.008 

-0.0125 

-0.0125 

-0.0125 

-0.013 

Slope / Falcon 

0.078 

0.052 

0.051 

0.0542 

0.0584 

0.063 

0.065 


Table 5. Values of the observed and simulated 0 3 /C0 slopes (JO3/JCO) under several 
conditions (with and without mixing, with and without clouds and with increased PAN 
decomposition rate. 




Figure 1. Flight tracks of the DC8 (18 July), the BAel46 (20 July) and the Falcon (23 
July) (black dashed lines). Lagrangian matches along the flights are represented in red 
together with FLEXTRA backward trajectories from these matches (blue for the Falcon, 
pink for the BAe-146 and violet for the DC8). 
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Figure 2. DC8 measured CO (black line) and jN02 (blue dashed line) together with 
simulated jN 02 along the flight with (pink dashed line) and without (pink solid line) 
aerosols. A low cloud layer was simulated at the blue rectangle location, and the BB 
plume at the red rectangle location. The largest jN02 variations are due to change in 
altitude of the aircraft (see flight track shown in Fig. 3). 
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Figure 3. Aerosol scattering ratio as measured by the DC8 on-board lidar on the 18 
July 2004. A strong signal is seen at the same location as high CO in-situ data (18h30, 
19hl5) including the Lagrangian match (18h39 - 19h02) 
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Figure 4. CiTTyCAT simulations initialized with mean and standard deviation mea- 
surements in the DC8 plume are represented respectively as continuous black lines and 
dotted color lines (blue for mean - std, pink for mean + std), respectively. Vertical red 
lines correspond to measurements made in the plume by the BAe-146 aircraft (20 July, 
202 Julian days), and the DLR Falcon aircraft (23 July), 205 Julian days). For the water 
vapor and temperature panel, the continuous line represents the simulated water vapor 
and the dotted line the simulated temperature. Red and green vertical lines represent the 
water vapor and temperature measurements, respectively. 



Figure 5. Change in PAN and O3 after 1 day of fixed simulation with different 
temperatures at the location and time of the DC8 plume sampling 
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Figure 6. CiTTyCAT simulations including mixing with a mixing timescale of 20 
days during the first 2 days and 5 days during the last 3 days. Same caption as Fig. 4. 
Background concentrations are given in Table 4. 



Figure 7. Observed O3 versus observed CO from the DC8 match (black crosses), the 
BAel46 match (red triangles) and the Falcon match (blue diamonds). 
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Figure 8. Measured O3 versus CO from the DC8 match (black crosses) used as 
initialization for the CiTTyCAT model, simulated O3 versus simulated CO after 2 days of 
simulation (red triangles) and after 5 days (bine diamonds). The lines represent fits to the 
measured correlations in the DC8, BAe-146 and Falcon matches, a) simulation without 
mixing, b) simulation with mixing, see text for details. 



Figure 9. 10 day backward trajectories from Jungfraujoch on the 24 July (pink), OHP 

on the 24 July (red), the Falcon BB plume sampling (green), and the MOZAIC BB plume 
sampling (blue) on the 23 July. Crosses correspond to the initial trajectory release time, 
stars to the BAe-146 mean match time (20 July, 12h30) and diamonds to the DC8 mean 
match time (18 July, 19h). 
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Figure 10. Measured CO (black line) and O 3 (pink line) at Jungfraujoch during July. 
Influence of BB plume are seen from the 24 July to the 26 July 


Figure 11. 
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Simulated 0 3 , CO, PAN, water vapor (continuous line) and temperature 


(dashed line) in the plume over 10 days with (pink) and without (blue) mixing. Mea- 
surements made by the BAel46 on the 20 July (red), the Falcon (red) and MOZAIC 
aircraft (green) on the 23 July, the OHP lidar (only 0 3 , orange) on the 24 July, and at 
Jungfraujoch on the 24 July (violet) are represented as vertical lines). 
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